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^ , We prove that spatial Kerr solitons, usually obtained in the frame of nonlinear Schrodinger equa- 



tion valid in the paraxial approximation, can be found in a generalized form as exact solutions of 
Maxwell's equations. In particular, they are shown to exist, both in the bright and dark version, as 
linearly polarized exactly integrable one-dimensional solitons, and to reduce to the standard paraxial 
form in the limit of small intensities. In the two-dimensional case, they are shown to exist as az- 
^ • imuthally polarized circularly symmetric dark solitons. Both one and two-dimensional dark solitons 

exhibit a characteristic signature in that their asymptotic intensity cannot exceed a threshold value 
C/3 , in correspondence of which their width reaches a minimum subwavelength value. 



I. INTRODUCTION 

The analytic description of spatial Kerr solitons, initiated by the seminal paper of Chiao et al.,P] has been con- 
tinuously evolving in the last forty years. It basically hinges upon the use of the nonlinear Schrodinger equation 
(NLS), which in turn follows from the nonlinear Helmholtz equation once the paraxial approximation, limiting the size 
a of the propagating beam to values large compared to the wavelength A, is introduced. This approximation becomes 
, inappropriate if the beam size a is comparable with A, a regime where nonparaxial effects become important and 
T-H ' are eventually able to provide a mechanism for avoiding nonphysical behaviors (like^e.g., catastrophic col lapse) i n the 
\ beam evolution. Although many contributions have been produced in this direction. |1. IsHo. ItL la. lol Hol fill flSj they 
' are typically based on some form of asymptotic expansion in the smallness parameter r) = X/a and are thus limited to 
c/3 I the range r/ < 1. To overcome this limitation, we start ab initio from Maxwell's equations and look for exact soliton 
. ^ ' solutions. More precisely, we solve Maxwell's equations in the presence of a fully vectorial Kerr polarizability and find 
a class of perfect optical solitons which inherently include all nonparaxial contributions. This is separately performed 
for one dimensional and two dimensional spatial solitons, both in the bright and dark configuration. In particular, 
the one dimensional case is dealt with by reducing Maxwell's equations to a system of first order differential equations 
and handling it by appealing to the usual formalism employed in the frame of dynamical systems. Our system is 
shown to posses a first integral so that its integrability is proved and the boundary value problem, associated with 
solitons, solved in closed analytical form. 

One the main results obtained in this paper is the proof of the existence of exact solutions of Maxwell equations 
d ' in the form of linearly polarized one dimensional Kerr solitons: they do not suffer of any limitation on the value of 
a and A (apart from the obvious ones associated with the validity of the macroscopic model of Kerr polarizability) 
and their existence curve can be numerically evaluated for all values of the beam intensity. Both bright and dark 
solitons can be derived from an integrable system of equations and their existence curve shows that, in the case of 
bright solitons, any value of the peak intensity u^q is allowed, while dark solitons can only exist if their asymptotic 
intensity does not exceed a threshold value completely determined by Kerr coefficients. In correspondence to 
this threshold, their width approaches the minimum value of the order of a fraction of A. In the two dimensional case, 
dark azimuthally polarized solitons are found and their existence curve implies, the same threshold behavior of the 
one dimensional dark solitons. While one dimensional solitons reduce to the standard paraxial ones for small values 
of the intensity, the two dimensional azimuthal dark soliton is a completely new entity which has never been studied 
in the paraxial regime. 
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We wish to note that the proof of the existence and derivation of exact sohtons requires, in the one dimensional 
case, the use of a rather sophisticated mathematical analysis borrowed from the dynamical system formalism, which 
we decide to report in full in section 2. 



II. ONE DIMENSIONAL SPATIAL SOLITONS 



The electric and magnetic complex amplitudes E(r) and B(r) of a monochromatic electromagnetic field 
-Re[Eexp(— iwi)], i?e[B exp(— iwt)] propagating in a nonlinear medium obey Maxwell's equations 



V X E = iwB, 



V X B 



(1) 



where rin labels the linear refractive index and Pni is the nonlinear polarizability. In the case of nonresonant isotropic 
media [lj|, the vectorial Kerr effect is described by the polarizability 



|E|2e+-(E-E)E* 



(2) 



n2 being the nonlinear refractive index coefficient. After eliminating B from Eq.JQ) and takind advantage of Eq.Q, 
we get 



V X V X E = k^E + P 

3 no 



lEpE- 



1 



(E • E)E* 



(3) 



where k = now/c. We now introduce a Cartesian reference frame Oxyz with unit vectors ex,ey,ez, and look for one 
dimensional solitons propagating along the z— axis, that is for y— independent fields of the form 



E(x, y, z) = exp(iQ;z) [Uxix)ex + iUz{x)ez 



(4) 



where Ux and Uz depends on x alone and a is a real constant. Substituting Eq.l^J into Eq.© yields the system of 
ordinary differential equations 



a— 

dx 

(fUz dUx 



dx^ 



dx 



no 



1, 



1 + + 
no V3 



Ux. 
Uz 



(5) 



whose unknowns Ux and Uz are real (as a consequence of the 7r/2 phase difference we introduced between the transverse 
and longitudinal field components (see Eq.Q). Note that the field in Eq.Q has a vanishing ?/— component, a 
requirement not forbidden by Maxwell's equations. From Eqs.(|SJl, it is also evident that the z— component Uz 
only vanishes if Ux = ±1/ (?T-o/2n2)(a^/fc2 — 1) which describes a family of solitary plane waves rather than solitons 
[l5| . The fact that Uz does not generally vanish is a consequence of the vectorial coupling between transverse and 
longitudinal components which cannot be rigorously neglected when describing spatially nonuniform fields, like for 
examples solitons (From a physical point of view, this follows from the first Maxwell equation setting the divergence 
of the electric field). Note that the longitudinal component is usually neglected in the paraxial regime thanks to the 
slow variation of the transverse component as compared to the wavelength A = 27r/fc, a circumstance which allows to 
treat it as a perturbation for slightly nonparaxial beams 0|. In the present paper, we deal on equal foot with 
both transverse and longitudinal components and it is their simultaneous non- vanishing and coupling which allows us 
to find exact solitons. 

Equations 10 can be recast in a more symmetric form by differentiating the first one and consequently eliminating 
d^Uz/dx"^ (together with dUz/dx) from the second one, thus getting 



dUx ^ [f3^ (1 - §7"^ + 2jul) + l{j + 2ul + lul) uj] 
/3[l+7(6«^ + H)] 



QxiUx,Uz\f 
Qz{Ux,Uz\, 



(6) 
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where we have introduced the dimensionless variables ^ = kx, f3 = a/k and {ux,Uz) = ^/l^iTl/^Ct^a;, t^z), while 
7 = "-2/1^12! (so that 7 = +1 and 7 = —1 for focusing and defocusing media, respectively). Equations © are a system 
of first order differential equations describing any electromagnetic field of the form of Eq.(@J and they are equivalent 
to Maxwell's equation, providing the relation 

1+7(6^2 + ^2^^0 (7) 

uniformly (i.e. for any ^) holds. Equations (|SJ) can be conveniently regarded as an autonomous dynamical system 
(since Qx and Qz does not explicitly depend upon whose solutions, or orbits, are ^—parameterized curves u(^) — 
{uxiO'^ziO)^ (belonging to the phase plane {ux,Uz)), tangent at each point to the vector field Q = {QxQz)^ ■ Solitons 
are particular orbits which, for suitable values of /3, pass through two special points of the phase plane imposed by 
the boundary conditions pertinent to each soliton kind (boundary value problem). 

The most remarkable and general property of the system of Eqs.@ is that it is a conservative system, i.e., it admits 
a first integral F[uxtUz\I3)^ defined over the phase plane, satisfying the relation 

^_dF _ dF dux ^ dF duz _ g + Q (8) 
dux d^ duz d^ dux duz 

In fact, it is straightforward to prove that the function 

F{ux.uz\p) = 2ul + ^utul + \ulu\ - 17(3/32 - 4)4 + ^7(2 - ti')ulul + \^0'u\ - + \l5'ul (9) 

obeys Eq.© whenever Eq.© is satisfied. This implies that F is a first integral of the system of Eqs.® whenever 
this system is equivalent to Maxwell's equations. According to Eq.®, any solution of Eqs.® is constrained to move 
along a single level set 

F(ux,Uz\P) = F^. (10) 

Inverting Ea. (|10|) furnishes Uz — Uz{ux, Fq, f3) which, once inserted into the first of Eqs.®, yields a first order 
differential equation solvable by quadratures, thus proving the integrahility of Eqs.®. Note that the first integral in 
Eq. (PJ) contains even powers of Ux and Uz only so that any level set of Eq. I|10|l is invariant under the inversion of the 
phase plane {ux,Uz) -{ux,Uz). 

Exploiting the properties of the first integral found above, we are now in a position to solve in a direct way the 
soliton boundary value problem, that is to find suitable values of f3 (if any) for which a solution Ux{£,)i Uz{£,) of Eqs.® 
satisfies the general boundary conditions 

Ux{0) \ _ f "2:0 \ _ „ f Ux{+oo) \ _ f 

^XOO 
" ''ZOO 



uzio) ; - V / - I uz{+oo) ] - Ur ' = ^^^^ 



where Uq and Uqo are defined by the kind of soliton, bright or dark, we wish to consider. From a geometrical point 
of view, this implies that the associated integral curve on the phase plane (ux,Uz) has to pass through the points Uo 
and Uoo, or, using Ea.({Tn|). 

F{UxO,Uzo\(3) = Fq, 
F{Uxoo,Uzoo\P) ^ Fq. (12) 

Since xioo has to be reached for ^ +00, it is obvious that xioo has to be an equilibrium point of Eqs.(|SJl, that is 

Qx{Uxoo,UzoG\li) = 0, 

QziUxca,Uzoc\l3) = 0. (13) 

Equations (|12|l and H13|l in the unknowns (3, Fq are necessary conditions for the solitons existence. They become also 
sufhcient if, once f3 and Fq are determined, one is able to prove that the integral curve actually reaches the point u^c- 
Following the outlined procedure, the existence of both bright and dark solitons will be proved and the corresponding 
existence conditions and propagation constants /? will be found. 
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A. Bright Solitons 

Bright solitons are localized nondifFracting beams, that is solutions of Eqs.lO vanishing for |^| —>■ +00, which in 
turn requires u^o = 0. Note that Eas. H13|l are automatically satisfied by this boundary condition since the origin 
{ux,Uz) = (0,0) is always an equilibrium point of Eqs.(|SJ|. The second of Eas. (|12|) directly gives Fq = so that the 
remaining condition we have to impose is the first of Eq. 112() that is 

F{u,n,u,n\(3)^0. (14) 

In order to set the boundary condition Uq we note that, because of the invariance of the level set in Ea. ()10|l under 
inversion of the phase plane, a soliton has to be associated with an integral curve starting from and ending into the 
origin and that this curve has to be symmetric under either the reflection — *■ —u^ or the reflection Uz — > —Uz- 
Because of these symmetry properties, we have uq = (0 Uzo)'^ and uq = (u^o 0)^ (where the symbol T stands for 
transposed) in the former and in the latter case, respectively. In the first case, Ea. ltTHl becomes /3^W2q(7u^o + 1) = 
which implies f3 — 0, so that soliton propagation is not allowed. We are left to consider the case Uo = (u^o 0)^ for 
which Ea. (|14l) furnishes 

1 + 3ju^ox 

In appendix A, we prove that bright solitons exist for all the real values of Uxo in focusing media (7 = 1) and that they 
never exist in defocusing media (7 = — 1) in agreement with the intuitive behavior of Kerr nonlinearity which tends 
to tighten and to spread the beam in these two cases, respectively. Obviously, the above results about the existence 
of bright solitons are based on the validity of Eq.©, which fails either for large intensities or for soliton widths so 
small to invalidate the continuum description of the material response. For 7=1, Eq. (|15|) yields 

^^^1 + 2^ (16) 
V 1 + 3uL 

which is the propagation constant of the exact bright solitons. The double sign in Ea. p6|l describes the two counter- 
propagating solitons along the z— axis. 

Substituting Ea. l|16|l and Fq = into Ea. Hl(J|) . we obtain the equation for the integral curves on the phase plane 
corresponding to bright solitons, and these are reported, for some values of Uxo, hi Figure 1. Note that, for each 
\uxo\j the corresponding level set is a bow-tie shaped curve encompassing three orbits of Eqs.®, that is the origin 
(which is an equilibrium point) and the left and right loop of the bow-tie. These last two orbits correspond to a pair 
of bright solitons each of which can be obtained from the other after the inversion of the x— axis, ^ — > — ^ (implying 
the reflection —Ux also), as expected because of the reflection invariance along any directions shown by Kerr 

nonlinearity. Considering the right half plane u^, > only, we observe that soliton curve Ux{S,), Uz{S,) explore the loop 
starting form the origin (for ^ = —00), reaching the point [ux^Uz) = (wa;o,0) (for ^ = 0) and ending into the origin 
(for ^ — -l-oo). From Eqs.® it is evident that the loop is explored counter-clockwise and clockwise for /3 > and 

/3 < 0, respectively, so that, for counter-propagating solitons (denoted with (-|-) and (— )), we have«i+^(0 =«i"'(0 
and ui+^(0 = 

Having proved the bright soliton existence and derived the associated propagation constant /3, we are now in the 
position to obtain the soliton shape for any given u^o by numerically solving Eqs.l© with [3 given by Ea. Hlti|l and the 
initial conditions Ux{0) = Uxq, Uz{0) = (the numerical approach being much simpler than integrating the system 
Eqs.© by quadrature). In Figure 2, we report the plots of the transverse Ux and longitudinal Uz components of the 
bright solitons for the same Uxo as in Figure 1. Note that, as expected, the soliton width decreases for increasing 
Uxo, while the longitudinal component Uz increases. In Figure 3, we report the bright soliton existence curve, relating 
the FWHM {Alright) to lu^^ol- As \uxa\ decreases the width indefinitely increases and diverges for \uxo\ = 0; on the 
contrary, as \uxo\ increases, the width decreases monothonically approaching zero. 



B. Dark Solitons 



In the scalar approximation, dark solitons are nondiffrating beams vanishing at ^ = and approaching an asymptotic 
amplitude value for |^| —f +00. In our vectorial case, the natural extension of the previous definition is identified with 
soliton solutions with Uq — (0 Uzo)^ , and Uqo — {uzoo 0)^ (see Eas. (|ll|) '). In fact, the above boundary conditions will 
be proved to describe an exact dark soliton which, in the paraxial limit, reduces to the standard scalar dark one. 
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The chosen values of Uoc identically satisfy the first of Ea. (|13|l . The second of Eas. (|13|) implies, with the help of 
the second of Eqs.®, 

f3' = l + 2^ul^. (17) 
Substituting this value of together with the boundary conditions into Eas. H12|l . we get 



,4 _ „,„,2 „,4 



7^20 - (18) 



The first of these equations furnishes the value Fq of the first integral along the dark soliton integral curve. The 
second one is a necessary condition for soliton existence from which we immediately obtain 7 = —1, in agreement 
with the intuitive property that only dcfocusing media can support dark solitons. In Appendix B, we prove that dark 
< 1/6 only and that 

fT 

UzO = ± 



Fo = i (1 - 2ulJ ut^, (19) 



2; 00 

so that each soliton is completely specified by the value u^oo only. 

As in the case of bright solitons, the integral curves in the phase plane associated to dark solitons are given by 
Ea. (|10|l . with f3 and Fq given in Eas. H19(l . some of them being reported in Figure 4. For each \uxoo\ the level set is a 
closed curve encompassing four orbits of Eqs.® that is the two equilibrium points {—u^^, 0) and {u.^.^, 0) together 
with the two curves joining these two points in the upper and lower half plane, respectively. These last two orbits are 
associated to a pair of dark solitons having opposite longitudinal components. Limiting our attention to the upper 
half plane Uz > 0, the dark soliton curve U2;(^),U2(C) starts, for (3 > 0, from the point (— WxocO) at ^ = ~oo, reaches 
the point (0,Uzo) at ^ = and finally ends into the point {uxoo,0) at ^ = +00 (for /? < it is sufficient to invert 

For any given value of Uxoo (in the range \uxoo\ < l/\/6), the shape of dark solitons can be obtained by numerically 
integrating Eqs.® with /? given by the first of Eas. l|19|) and initial conditions Ux{0) = and Uz(0) = Uzo (the latter 
being given by the second of Eas. (|19l) ). In Figure 5, we plot the transverse Ux and longitudinal Uz components of 
various dark solitons, for the same Uxoc as in Figure 4. Also in this case, for increasing Uxoo the soliton width decreases 
while the longitudinal component increases. In Figure 6 we report the dark soliton existence curve relating the soliton 
FWHM (Adark) to Uxoo: the range < Uxoo < l/'v/S- Note that, for very small Uxoo- the FWHM indefinitely grows 
whereas in correspondence to the threshold value Uxoo = l/\/6, it attains its minimum value ~ 4, corresponding to 
dimensional value ~ (2/7r)A ~ 0.63 A. 



C. The Optical Intensity 



Having derived the electric field (see Eq.Q) associated to both bright and dark solitons, we can directly evaluate 
the corresponding magnetic field by means of the pertinent Maxwell equation. Substituting Eq.Q) into the first of 
Eqs.0J we easily deduce, in terms of the dimensionless fields. 



B(x, y, z) = - . / exp(z/3fcz) ( - ^ 
LU V \n2\ V 



^ — kx 



(20) 



Note that the soliton magnetic field is parallel to the y— axis and therefore orthogonal to the electric field everywhere, 
a remarkable vectorial feature that exact solitons shares with plane waves (which are rigorously nondiffracting fields 
as well). In order to describe the soliton energy flow, we can now evaluate the averaged Foynting vector S = 
i?e(E X B*)/(2/io) which, using Eqs.0) and if^ and the second of Eqs.®, becomes 



1 + 27 + 



1 



uIbz 



IGz, 



(21) 



where Iq = kno/ (2ujfio\n2\) and /, the modulus of the averaged Foynting vector, is the optical intensity. The averaged 
Foynting vector lies along the z— axis everywhere and this is fully consistent with the nondiffracting nature of the 
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solitons we are considering (which is not rigorously the case in the paraxial approximation) . Note that S is proportional 
to /3~^ and the expression in square brackets of Ea. (|21|l is always positive (while this is trivial in the case 7 = +1, 
in the case 7 = —1 all the orbits Ux{Ot''^z{£,) of Eqs.© lie inside the ellipse defined in Eqs.Q, which is in turn 
contained within the ellipse 2 (it^ + ^u^) = 1, so that the expression in square bracket of Ea. H21l) is always positive). 
This implies, as expected, that, both for bright and dark solitons, the sign of (3 determines whether S is parallel 
or antiparallel to the z— axis. In Figure 7 we report the plots of the normalized optical intensity I /Iq for the same 
bright and dark solitons examined in the previous Figures. From Ea. H21(l we observe that the optical intensity is in 
general not proportional to the square modulus of the electric field. However, in the paraxial limit where Ux << 1, 
Uz «Ux and ~ 1, Eq.|(2U gives / ~ Iqu^, reproducing the well-known result typical of paraxial optics. We can 
also evaluate the maximum soliton optical intensity, that is Ea. H21|) at ^ = (and 7 = +1) for bright solitons and at 
1^1 = +00 (and 7 = —1) for dark solitons, thus getting 

hright = ^0\/l + 3u^o"xO' 

Idark = lo Vl - 2u2^U^^ . (22) 

From these equations we note that Ibright > IqUxq whereas Idark < Iqu'^oo that, in general, bright and dark solitons 
are characterized by an optical intensity which is greater and smaller, respectively, than the corresponding paraxial 
prediction. This is evidently associated with the fact that, in an extremely narrow soliton, the longitudinal component 
of the electric field is as large as the transverse one. 



D. Paraxial limit 



The above description of one dimensional bright and dark solitons is exact, no approximation having been exploited 
in their analytical derivation. As a consequence, the solitons described above are expected to reduce, in the paraxial 
limit where the soliton width is much larger than the wavelength, to those predicted by the NLS. The paraxial limit 
is clearly obtained by considering the range of values 

Ux « 1, 

Uz « Ux, (23) 

since the soliton width increases for decreasing optical intensities while the longitudinal component decreases. By 
differentiating the first of Eqs.®, using the second of Eqs.® to eliminate duz/dS^ and exploiting Eas. (P5|l to retain 
only the first relevant order, we obtain 

^ = - l)u. - 2^ul. (24) 

Note that, in describing paraxial Kerr solitons, the electric field is usually expressed as Ex{x,z) — exp[zA:(l + 
f3)z]^/noJ\n2\ux{£,), where the fundamental plane wave carrier exp(ifcz) is separated by the slowly varying ampli- 
tude of the field. The comparison of this field expression with Eq.Q) yields /3 = 1 + P with (3 << 1 (consisting with 
the paraxial picture where the main plane wave carrier is slowly modulated), implying that — 1 ~ 2/3. Introducing 
this relation into Ea. H24|l . we get 

~ 1 d'^Ux Q 

2^ = (25) 

which coincides with the usual equation (obtained from the NLS) describing paraxial Kerr solitons. Equation H25|) 
admits both of bright soliton solutions (7 = +1) of the form Ux{£,) = Uxo sech(ua;o^) and of dark soliton solutions 
(7 = —1) of the form Ux{^) — Uxoo tanh(wa;ooO- The propagation constants are respectively given by /? = u'^q/2 and 
j3 = —UxQOJ which can also be found, mutatis mutandis, from Eqs.(16) and the first of Eqs.(19), whenever the paraxial 
conditions {uxo << 1 for bright solitons and Uxoo << 1 for dark solitons) are satisfied. These solitons obviously 
coincide with the asymptotic paraxial limit of the solitons described in this paper. In order to make this comparison 
more quantitative, in Figure 3 and Figure 6 we have superimposed to the exact soliton existence curves (solid curves) 
their paraxial counterparts (dashed curves). More precisely, the FWHM of bright and dark paraxial solitons are 
easily shown to be Alright — 2.6348/Ma;o and A^ark — 1.0986/ua;oo- As expected, the paraxial and exact curves are 
practically indistinguishable for small Uxo or Uxoo- Not surprisingly, for dark solitons, the agreement between exact 
and paraxial prediction is satisfactory almost everywhere since the value of Uxoo is restricted to be less than l/\/6 
that is to a moderate nonparaxial regime. 



7 



III. TWO DIMENSIONAL CASE: AZIMUTHALLY POLARIZED SPATIAL DARK SOLITONS 



In order to deal with the two dimensional case, we introduce polar cylindrical coordinates r, (p, z with unit vectors 



er^e^p^ez and look for fields of the form 



E(r, ipz) = E^{r, z)e^ + Ez{r, z)ez, 



(26) 



describing a circularly symmetric configuration with vanishing radial component. Inserting Eq.ljJHJ in Eq.Q, we 
obtain 



drdz 



0, 



d'^E.^ , d fdE, 
dz'^ 
d'Ez 



dr \ dr r 



= -eE^-e^^' 



3 no 



:(E-E)K 



IdE, 



dr 



= -k'Ez 



3 no 



m^Ez + -(E--E)E: 



(27) 



Internal consistency of the set of Eas. l|27|l (three equations in two unknowns) requires E^ = 0. As a consequence, the 
second of Eas. (|?7| yields 



d^^d_ fdE^ ^ E^ 
dz"^ dr \ dr r 



-k^E^^2k^ — \E^\^E^. 
no 



(28) 



We note that circular symmetry and polarization imposed to the field, together with the symmetry properties of Kerr 
effect, have allowed us to reduce Maxwell's equations to the single Ea. H28|) . Equation (|28|l is conveniently rewritten 
in the dimensionless form 



dC^ dp \dp p 



= -U-2j\U\^U, 



(29) 



where p = ^/2kr, ^ — kz, U = \/\n^\/noEip and 7 = n2/\n2\. If we look for soliton solutions of the form 

Uip,0^e^''<uip), (30) 

Eq.(|25J) becomes 



— ( — + - 1 = -(a^ - l)u - ^u^ 
dp\dp p' ' ' 



(31) 



Both the structure of Ea. (|31|) and the azimuthal field polarization dictate u(0) — 0, so that azimuthally polarized 
bright solitons do not exist. In order to find dark solitons, we introduce the further condition 



lim u{p) = Uoo, 

p — ^00 



(32) 



together with the vanishing of all derivatives for p —^ 00. Since focusing media (7 = 1, i.e., n2 > 0) are not able to 
support dark solitons, we consider hereafter defocusing media (7 = —1, i.e., ^2 < 0), so that Ea. (|^ reads 



d f du u\ 1,9 , Q 
which implies, together with the above boundary condition at infinity. 



(33) 



(34) 



While positive and negative signs of a respectively refer to forward and backward travelling solitons (see Ea. (|3(Jf) '). 
u[p) depends on (see Ea. (|^ '). Equation shows the existence of an upper threshold for the soliton asymptotic 
amplitude 



V2' 



(35) 
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since, otherwise, a would become imaginary. If we now insert Ea. H34(l into Ea. (|33() . we obtain 

We have carried out a numerical integration of Eg. 1)36(1 with boundary conditions u(0) = and u{oo) = Uoo, by 
employing a standard shooting-relaxation method for boundary value problems. The results of our simulations show 
that dark solitons can be obtained in the range of field amplitudes < Uqo < 1/ \/2- Different soliton profiles are 
reported in Figure 8. 

In order to complete our analysis, we now evaluate both the magnetic field and the Poynting vector. Recalling the 
expression of the soliton electric field 



E = J^e'"'''u{V2kr)e^, (37) 

1*^2 I 



we obtain, from the first of Eqs.(^ written in cylindrical coordinates 

712 W 



aufir + iV2 , 

dp p 



(38) 



The magnetic field has a radial component whose shape coincides with that of the electric field, and a vanishing 
azimuthal component, so that E and B are mutually orthogonal. With the help of Eas. (|37|l and 138() . the time 
averaged Poynting vector 

S = -^i?e(E X B*) (39) 
2po 

turns out to be given by 

Sir) = ^^.2(V2fcr)e. = ^|Epe.. (40) 
2ujpo \n2\ 2Ljpo 

We note that S is parallel to the z— axis, consistently with the shape-invariant nature of solitons. From an analytical 
point of view, this corresponds to the 7r/2 phase difference between and E^p (see Eas. (|37|) and (|38l) ). As expected, 
the Poynting vector is either parallel or antiparallel to according to the sign of a, while its amplitude is proportional 
to |Ep. The above plane wave-like properties are consistent with the nondiffractive nature of exact solitons. 

It is worthwhile to underline that, in the case of the azimuthal dark solitons we are considering, the asymptotic 
optical intensity loo — |S(oo)| turns out not to be proportional to u^. In fact, by using Eas. (|^ and one obtains 



/oo(Uoo)=/0"LVl^2^ (41) 

where Iq ~ fcno/(2a;/io|n2|), whose profile is reported in Figure 9. Equation (|4f |l shows that the asymptotic optical 
intensity is not a monotonically increasing function of the asymptotic field amplitude, but reaches its maximum 
threshold value I™""^ = /q/S'^/^ in correspondence to Mqo = l/\/3- This is connected to the a— dependence of the 
magnetic field (see Ea. ((38|l ) whose radial part tends to vanish for Uoo i/V2- A related and relevant consequence 
of Eq. H41|) is the existence of two solitons of different widths for a given asymptotic optical intensity. The existence 
curve relating the normalized half width at half maximum (HWHM) of the soliton optical intensity profile |S(p)| to 
Moo is reported in Fig. 10. In particular. Fig. 10 shows the existence of a normalized minimum HWHM ~ 2.1 (~ 0.24A) 
for Moo ~ l/\/2- In addition. Fig. 9 shows that a normalized HWHM ~ 2.7 (~ 0.3A) corresponds to lioo = l/\/3 for 
which the soliton attains the maximum asymptotic optical intensity /™"^. 

It is interesting to examine the behavior of our solution in the limit of large p. To this end, neglecting in Ea. l|36|) 
the term in u/ p, we have 

= (-^-"D" (42) 

which formally coincides with the equation describing one dimensional linearly polarized paraxial dark solitons. Equa- 
tion (|42l) admits of the solution u = Uoo tanh(/5Uoo/'\/2)- This solution can be compared with the exact one. This is 
done in Figure 11 where the ratio R{p) between the hyperbolic tangent and the exact solution is reported as function 
of p, for different values of Moo- The hyperbolic tangent solution reproduces the exact one for large values of p, as 
expected, while it at most differs by a factor = 1.2 for small values of p. 
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IV. CONCLUSIONS 

In this paper, the problem of the existence of nonparaxial spatial Kerr solitons has been completely solved. This has 
been done by showing that spatial solitons can be derived as exact solutions of Maxwell equations (thus making, within 
our approach, the term "paraxial" and "nonparaxial" redundant). In the one dimensional case, the perfect optical 
soliton represents the straightforward generalization of the paraxial one, the main difference being that dark solitons 
exhibit, unlike their paraxial counterpart, a specific upper limit for the possible values their asymptotic intensity can 
assume. In the two dimensional case, the exact dark soliton is of a completely new kind, and the difference between 
paraxial and nonparaxial becomes rather meaningless. In any case, the comparison between paraxial and exact 
solitons, done, for example, by inspecting the relative existence curves, shows that our solitons are a definite entity, 
independent from used the approximation scheme: the transition between the paraxial and the highly diffractive 
regime is very smooth and does not exhibit any kind of dramatic catastrophic behavior, as implied by the standard 
paraxial theory. 

APPENDIX A: APPENDIX A: EXISTENCE OF BRIGHT SOLITONS IN FOCUSING MEDIA 

In order to tackle the problem of bright solitons existence we have to prove that the curve defined in Ea. (|14|l . with 
given by Ea. (|15|l . actually reaches the origin of the phase plane {ux^ Uz)- To this end, it is convenient to introduce 
the polar coordinate defined hy = p cos (j) and Uz — p sin 0, so that Eq. H14|l becomes 



+7P 



p^ \ p^ COS^ ( 



2C0S'' 



4 9 ( ■ 2 y ^ • 4 

' -h - COS <p sm -h 77 sm i 



^(4 - 3/3^) cos'^ </> + i(2 - cos^ <?!)sin^ <i> + \^0^ sin"* , 



+ i [(1 - i?) cos^ + sin^ I = 0. 



(Al) 



This equation is trivially satisfied by setting p — ^ and this is consistent with the fact that the origin is in itself an 
orbit of Eqs.(|SJ|. Therefore, the integral curve associated with solitons is described by the vanishing of the expression 
within the curly brackets. Requiring that this curve reaches the origin yields 



tan 



1 



1 



'■ax 



4Mfa + 7 
(1 + 27^2 J 



(A2) 



where 0o = 4'{p = 0) and 0^ has been obtained from Eg. ljlSII . For 7=1, the RHS of Ea. ljA2p is positive so that 
this equation can always be solved which, together with the fact that Eq.© is always satisfied for 7 = 1, implies 
that bright soliton exist in focusing media for any value of uoa:- In the case 7 = — 1, the RHS of Ea. (jA2|) is positive 
for luorcl > 1/2 so that, the curve actually reaches the origin. However, in this case, the curve joining the points 
{ux,Uz) — (0,0) and {ux,Uz) = {uxn,Q) unavoidably crosses the ellipse 6u^ + (2/3)m^ — 1 since its semi-axis along 
the X— axis is l/\/6 < 1/2. Therefore, for 7 = —1, a point belonging to the integral curve such that Eq.Q fails to be 
satisfied always exists, with the consequence that, in defocusing media, bright solitons never exist. 



APPENDIX B: APPENDIX B: CONDITIONS FOR DARK SOLITONS EXISTENCE 

As already explained in Section 2, Eas. ljl?!) and (|18|l are necessary for dark soliton existence so that we have to 
find when they are also sufficient. From Ea. (|17|) (with 7 = —1) it is evident that solitons can exist for < 1/2. 
The equation for the dark soliton integral curve on the phase plane fEa. HlU|) with f3 and Fq given in Eas. H19|) ) can be 
solved for thus yielding 



2 32; 
= 



9 "2: 



2(1 



9 "a: 



4(1 



■ 2u2^)] [4u6 - (1 + 6ul^)ui + 2ul^ul - (1 - 2ul^)u\ 



" 1 

2;ooJ 



9 "a: 



2(1 - 2^2, 



(Bl) 
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which furnishes as a function of (parametrically dependent on u^^) along the dark sohton integral curve. 
Here, the plus sign between the two terms has been chosen in order to satisfy the boundary condition u^iu^^) = 0. 
Evaluating Ea. (|Bl(l at = and taking the square root of the result we obtain the second of Eas. H19|l . which is 
consistent with the boundary conditions since it satisfies the second of Eas. (|18|l (with 7 = —1). Therefore, in order to 
prove soliton existence, we are left with proving that the curves in Ea. (|Bl|) actually reach the point Uqo — {uxoo 0)"^ 
(i.e., with proving that the RHS of Ea. l|Bip is a positive real number) and that such curves do not cross the ellipse 
+ {2/3)u1 = 1, thus leaving Eq.© satisfied. Since the expression under the square root is always positive for 
■"^oo < 1/2, we have only to ensure that the RHS of Ea.(|Bl|l is positive. It is not difficult to show that this is the 
case whenever 

Aul - (1 + 6ul^)ul + 2ul^ul - (1 - 2ul^)ul^ < 0. (B2) 

Imposing that the maximum of the polynomial in the LHS of this inequality is negative, we obtain the condition 
^xoo < 1/*^- The existence of dark solitons in this range for u^^o finally proved by noting that any integral curve 
associated to these solitons globally lies within the ellipse + (2/3)m^ = 1, so that Eq.((7| is always satisfied. 
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Figure Captions 

• Figure 1: Plot of phase portrait of Eqs.© associated to bright solitons for \uxo\ — 1,2,3,4. Each bow-tie 
shaped curve is obtained by plotting the level set defined in Ea. fT^ with f3 given by Ea. lfTHji . Any piece of 
curve starting from and ending into the origin (left or right loop of each bow-tie) is associated to a single bright 
soliton. 

• Figure 2: Plot of the transverse component Ux{£,) (a) and longitudinal component Uz(^) (b) of bright solitons 
for Uxo = 1, 2, 3, 4 (same cases as in Figure 1) and f3 > 0. 

• Figure 3: Bright soliton existence curve (solid line), relating the FWHM, Alright, of the amplitude Ux{^) to 
|u2,o|- For very small and very large \uxo\, the FWHM diverges and vanishes, respectively. The dashed line 
represents the FWHM, Alright, of paraxial bright solitons. Note the complete overlapping of the two curves for 
Uxo < 0.2. 

• Figure 4: Plot of phase portrait of Eqs.© associated to dark solitons for \uxoo\ = 0.1,0.2,0.3,0.4. Each loop 
is obtained by plotting the level set defined in Ea. H10() with /3 and Fq given in Eas. (|19|l . Any piece of curve 
joining the points ( — Uxoo^ 

0) and {uxoo,0) is associated to a single dark soliton. 

• Figure 5: Plot of the transverse component Ux{(,) (a) and longitudinal component Mz(^) (b) of dark soltions 
for Uxoo = 0.1, 0.2, 0.3, 0.4 (same cases as in Figure 4) and /3 > 0. 

• Figure 6: Dark soliton existence curve (solid line), relating the FWHM, Adark of the amplitude Ux{^) to Im^ooI- 
For very small lu^^ool, the FWHM diverges whereas at the threshold value Uxoc = l/'v/S it attains its minimum 
value ~ 3. The dashed line represents the FWHM, A^ark, of paraxial dark solitons. Note the complete overlap 
for most of the values of Uxoo- 

• Figure 7: Normalized optical intensity |S(^)|//o of bright (a) and dark (b) solitons evaluated from Ea. l21() for 
the same soliton conditions as in Figure 1 (for bright solitons) and Figure 5 (for dark solitons). 

• Figure 8: Two dimensional dark soliton profile u{p) for various values of Uqo. 

• Figure 9: Normalized asymptotic optical intensity Ioc/Iq as a function of the asymptotic dimensionless field 
amplitude Uoc- Note that two solitons exist for any allowed asymptotic optical intensity. 

• Figure 10: Existence curve relating the normalized soliton optical intensity HWHM to Uoo- 

• Figure 11: Plot of the ratio R{p) — Uoota,nh{uoo p/V^)/u{p) for different values of Moo- 
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